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ABSTRACT 

The gravitational waves (GW) emitted during the coalescence of supermassive black holes (SMBHs) in the 
mass range ~ (1O 4 -1O 7 )M0/(1+z) will be detectable out to high redshifts with the future Laser Interferometric 
Space Antenna (LISA). The distance and direction to these "standard sirens" can be inferred directly from the 
GW signal, with a precision that depends on the masses, spins and geometry of the merging system. In a 
given cosmology, the L/SA-measured luminosity distance translates into a redshift shell. We calculate the size 
and shape of the corresponding three-dimensional error volume in which an electromagnetic counterpart to a 
LISA event could be found, taking into account errors in the background cosmology (as expected by the time 
LISA flies), weak gravitational lensing (de-)magnification due to inhomogeneities along the line of sight, and 
potential source peculiar velocities. Weak lensing errors largely exceed other sources of uncertainties (by a 
factor of ~ 7 for typical sources at z = 1). Under the plausible assumption that SMBH-SMBH mergers are 
accompanied by gas accretion leading to Eddington-limited quasar activity, we then compute the number of 
quasars that would be found in a typical three-dimensional LISA error volume, as a function of BH mass 
and event redshift. Low redshifts offer the best opportunities to identify quasar counterparts to cosmological 
standard sirens. For mergers of ~ (4 x 10 5 - 10 7 )M Q SMBHs, the LISA error volume will typically contain a 
single near-Eddington quasar at z ~ 1. If SMBHs are spinning rapidly, the error volume is smaller and may 
contain a unique quasar out to redshift z ~ 3. This will allow a straightforward test of the hypothesis that 
GW events are accompanied by bright quasar activity and, if the hypothesis proves correct, will guarantee the 
identification of a unique quasar counterpart to a LISA event, with a B-band luminosity of Lb ~ (10 10 - 10 n )L©. 
Robust counterpart identifications would allow unprecedented tests of the physics of SMBH accretion, such as 
precision-measurements of the Eddington ratio. They would clarify the role of gas as a catalyst in SMBH 
coalescences, and would also offer an alternative method to constrain cosmological parameters. 
Subject headings: cosmology: theory - cosmology: observations - large scale structure of universe - cosmic 
microwave background- galaxies: clusters: general 



1. INTRODUCTION 

One of the main objectives of the Laser Interferometric 
Space Antenna (LI SA), to be launched around the year 2013 
( Danzmannl l2004l) . is to detect the gravitational wave (GW) 
signals associated with coalescing supermassive black holes 
(SMBH) at cosmological distances. The LISA detector is de- 
signed to be particularly sensitive in the frequency range be- 
tween (3 x 10" 5 - 10" 4 )Hz < / < 0.1Hz, allowing the detec- 
tion of binary coalescences with total masses between 10 4 and 
10 7 M Q out to high redshifts. The limiting redshift depends 
on several factors (such as the orientation of the spins and or- 
bital plane of the SMBH binary, and its location on the sky 
relative to LISA ), and is expected to lie between z ~ 5- 10 
(Hushes 2002). A comparison of the gravitational wave- 
form with the anticipated detector noise can be used to esti- 
mate the accuracy with which LISA will be ab le to extract the 
physical parameters of the coalescenc e events (|Hugh£s |2Q02|; 
iBarack & Cutlerl2004llVecchiol2004llHolz & Hughes I200l . 



Of particular interest, in the context of searching for electro- 
magnetic (EM) counterparts, is whether the spatial location of 
the GW event can be localized to within a sufficiently small 
three-dimensional volume. In this paper, we determine the 
probability of finding a unique EM counterpart within the ex- 
pected error volume associated with SMBH merger events, 
for a range of masses and redshif ts. 

It has been argued by Vecchio (2004) that the identification 
of such EM counterparts will be difficult because, in typical 
cases, there will be too many co unterpart candidates to choose 
from. However, Vecchio (2004) associated counterparts with 
host galaxies and galaxy clusters, and used only the 2D an- 
gular positioning information for the analysis. In contrast, 
here we account for the 3D spatial information by using the 
redshift of an electromagnetic counterpart candidate in rela- 
tion with the luminosity distance determined by LISA and we 
focus on quasars as plausible counterparts. With these speci- 
fications, we shall demonstrate that in some cases, a specific 
counterpart can be uniquely determined. 
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If electromagnetic counterparts to LISA events exist, they 
will likely be related to the accretion of gas onto the SMBHs 
involved in the coalescence. Provided this accretion is not 
supply-limited, bright quasar counterparts approaching the 
Eddington luminosity would then be expected. A few ad- 
ditional arguments favor this scenario: galaxy mergers in 
hierarchical scenarios of structure formation are expected 
to deliver a significant amount of gas to the central re - 
gions of the merging galaxies ( Barnes & Hern auistl 119921) . 
and this g as may play a catalyst role in driving SMBH co- 
alescence f Begelman. Bla ndford. & Reeslll980l:lGould & Rixl 
120001 lEscala et al j|2004). Ultimately, however, many of the 
complex processes involved remain po orly understood. For 
example, Armitage & Nataraian (2002) have argued that, in 
the limit of a small mass ratio of the two SMBHs, a prompt 
and luminous electromagnetic sign al may be expe cted during 
coalescence, while Milosavlievic & Phinnev (2005) have ar- 
gued that in the limit of equal mass SMBHs, only a much 
delayed electromagnetic afterglow would be expected. All 
"cosmological standard sirens" may thus not be equal in their 
potential for electromagnetic counterparts. 1 Our working as- 
sumption in the present study is that bright quasar activity is 
a plausible electromagnetic counterpart to LISA events. This 
allows us to quantify the feasibility of an unambiguous iden- 
tification of such a counterpart. As we shall see below, the 
search for the counterparts will allow a test of the assumption, 
as well. 

The secure identification of the EM counterpart to even a 
single GW event could be useful in different ways: (1) to im- 
prove our understanding of the SMBH accretion physics, (2) 
to clarify the role of gas as a catalyst in SMBH mergers and 
(3) to supply an independent constraint on the background 
cosmology. A joint GW - EM analysis could, in principle, 
determine the masses and orbital parameters of the SMBH bi- 
nary, and yield a precise measurement of the Eddington ratio, 
L/LEdd) which will supply a key param eter in studies of the 
evolut ion of the BH/quasar population ( Small & Blandford 
19921 lHaiman & Loeti 119981 iKauffmann & HaehneU I2000t 
Haiman & Menoul 120001 IWvithe & Loebl 120031) . This pa- 
rameter is currently p oorly known (constrained by i ndirect 
empirical correlations; Vestergaard 2004, Woo & Urrv 2002, 
Kaspi ~al.l2000l) . The values range from ss 0.1 to > 1, with 
indications that higher-z quasars may be closer to LEdd than 
the z ~ quasars. Likewise, a joint GW - EM analysis could, 
in principle, b e used to estimate cosmological parameters 
dSchutzll986l) .bv comparing the luminosity distance (which 
is a direct observable by GWs) with the redshift (as inferred 
from the spectrum of the counterpart - in this case, a quasar). 
This would serve as a complement to constraints from the 
luminosity distance to high-z type la Supernovae (SNe), but 
with different systematic errors, and with the potential of ex- 
tending to higher redshifts. New constraints, spanning the 
range < z < 2, would be particularly well-suited to probe 
the properties of dark energy, which is expected to become 
dynamically dominant within this cosmic epoch. 

The rest of this paper is organized as follows. In § [2] we 
summarize our method to estimate the angular and radial po- 
sitioning errors expected from LISA, for SMBHs with a range 
of masses at different redshifts. In § [5] we discuss the con- 
version of the luminosity distance, as determined by LISA 
from the GW signal alone, to the redshift of the source. In 

1 The name "standard siren s" was suggested by Sterl Phinney and Sean 
Carroll <HoTz & Hughes 12555) . 



TABLE 1 
LISA Measurement Errors 





SM/M 




6d L /d L 


5U 


best 


0.8 x 1(T 5 


2 x 10~ 5 


2 x 10~ 3 


0.01 deg 2 


typical 


2 x 1(T 5 


9 x 10~ 5 


4 x lO- 3 


0.3 deg 2 


worst 


0.8 x 10~ 3 


0.1 


2 x 10~ 2 


3 deg 2 



NOTE. — Assumed SMBH binary parameters: mi = 
m 2 = 10 6 M Q andz = 1. 



particular, we discuss the uncertainty in the resulting redshift 
estimate. In § 0] we discuss our estimates for the number of 
quasars that may be found in the 3D error volume provided by 
LISA, based on the luminosity function and clustering proper- 
ties of known optical quasars. In § |5j we present our main 
results, and show that for typical low-redshift GW events dis- 
covered by LISA, a unique quasar counterpart may be identi- 
fied. In §|6] we point out various implications of a successful 
identification and discuss several caveats to this conclusion. 
Finally, in §0we summarize our conclusions. Unless stated 
otherwise, throughout this paper we assume a standard cold- 
dark matter cosmology (ACDM), with (ft\, VLm, Hq) = 
(0.70, 0.30, 0.047, 70 km s" 1 Mpc" 1 ), consistent with the re- 
cent result s from the Wilkinso n Microwave Anisotropy Probe 
(WMAP){ Sperg el et alJl2003l) and the Sloan Digital Sky Sur- 
vey (SPSS') llTe gmark et al 120041) . 

2. LOCALIZING LISA EVENTS 

A few studies have been carried out so far to address how 
accurately LISA will measure the source parameters of a co- 
alescing pair of SMBHs. In general, the accuracy depends 
on a large number of parameters: a total of 17 parameters 
in the most general case include 2 red-shifted mass param- 
eters, 6 parameters related to the BH spin vectors, the or- 
bital eccentricity, the luminosity distance, 2 angles identify- 
ing the position on the sky, 3 angles that describe the orien- 
tation of the orbit, a reference time, and a reference phase. 
Pue to the resulting computational limitations, various stud- 
ies have concentrated on small portions of the parameter 
space. The most up-to-date calculations estimating param- 
et er un certainties for SMBH i n-spirals have been carried ou t 
bv lBerti Buonanno. & Wilil(l2005al> : lHolz & Hugheslll2005h : 
and by | Vecchiol 1120 04). As compared to previous stud- 
ies, IVecchiol yp04) accounts for the effects of spins, and 
shows that parameter estimation errors improve significantly 
(by a factor of 3-10 for high spins) for selected parame- 
ters. IVecchiol {2004) also adopts an optimistic LISA sensi- 
tivity curve, by adopting the smallest observable frequency to 
be ~ 3 times lower than previous estimates and only considers 
cases with equal mass SMBHs. Our analysis, which relies on 
Vecchio's estimates, is therefore approximate to this extent. 

For con creteness , we ad opt the parameter uncertainties ob- 
tained by Vecchio (2004) for an equal-mass SMBH binary 
with m\=mi = 1O 6 M at redshift z = 1 . The uncertainties vary 
as a function of the fiducial orientation of the source relative 
to LISA, and are primarily influenced by the BH spin magni- 
tudes (i.e. higher spins lead to smaller uncertainties). Here 
we distinguish three cases. For our "best" case, we adopt the 
errors that correspond to the 10th percentile of the distribu- 
tion of uncertainties obtained by Vecchio ( 2004) for high BH 
spins (with dimensionless spin a = S/M 2 = 0.9, where S is the 
magnitude of the total spin and M is the total mass). For our 
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"typical" case, we adopt the errors corresponding to the 50th 
percentile in the case of moderate BH spins, with a = 0.5. 
For our "worst" case, we adopt the 90th percentile of the no 
spin case (a = 0). Note that SMBHs are ge nerally expected to 
be spinning fast (e.g. Volonteri et al. 2005), so that our "best" 
case may actually be representative of a fair fraction of events. 
In Tabled we list the errors on the chirp and reduced masses, 
M. = (m\miy'l i /{m\ +m2) 1//5 and /i = (jnxmi) / ! {m\ + rri2) re- 
spectively, and on the GW source location (i.e. luminosity 
distance, c/ L , and solid angle, £1), for these three cases of in- 

teres t. 

In IVecchiol (12004) . parameter errors have only been esti- 
mated for a single choice of an equal-mass SMBH binary 
with total mass Mq = m\ +ni2 = 2 x 10 6 M Q and redshift zo = 1 ■ 
Starting from this result, we crudely estimate the uncertainties 
Wl and 5il for other combinations of masses and redshifts as 
follows. First, note that the luminosity distance is simply pro- 
portional to the inverse of the signal amplitude. Therefore its 
estimator depends primarily on the total signal power, rather 
than on the specific shape of the signal waveform 2 . The lu- 
minosity distance error would then obey the simple scaling 



6d L (M,z) 



SNR(M lZ ) 
SNR(M ,z ) 



6d L (M ,zo) 
d L (M Q ,z ) 



(1) 



where SNR(M,z) is the expected value of the signal to noise 
ratio of the detection, 



•// a (*f t ,Ar) JiAj) 



(2) 



Here, M z = (1 +z)M is the red-shifted total mass, h(f,M z ,z) 
denotes the Fourier decomposition of the signal detected 
by LISA, and S„(f) is the RMS noise density per fre- 
quency interval (includ ing instrumental and confusion noise, 
iBarack & Cutlerlf2004P '). A crucial parameter for high M z 
values inscribed in S„(f) is the low frequency wall of the 
detector, which is further discussed below. The integration 
bound /i SC0 (M z ) corresponds to the innermost stable circular 
orbit (ISCO), beyond which the gravitational waveform is not 
well known, and f a (M z ,AT) is the arrival frequency a time 
AT before the ISCO is reached by the coalescing binary (see 
Vecchio 2004 for details). Throughout this paper, we fix the 
observation time of SMBH binaries at AT = 1 yr, unless a 
binary is so massive or at such a high redshift that it is not 
observable by LISA for a full year and then AT < 1 yr. Note 
that h(f,M z ,z) depends on other parameters, such as the an- 
gular momentum vector orient ation relative to LISA'S arms, 
the magnitude of spins, etc. (see Vecchio 2004 and references 
therein). For this estimate we calculate the SNR with the lead- 
ing order (i.e. Newtonian) contribution. The resulting depen- 
dence of the signal-to-noise ratio on BH masses and redshift 
are shown in Figure ^ The figure shows that the sensitiv- 
ity degrades significantly for distant sources, and also that it 
peaks in the mass range of 10 s - 10 6 Mq, which produce GWs 
near the optimal frequencies within the LISA band. 

The other important parameter, the angular position, is ex- 
tracted from the change in the relative orientation of LISA dur- 
ing its orbit around the Sun. Since LISA'S orbital time-scale 

2 The signal power also scales with the red-shifted chirp mass as A4^ 6 . 
However, this parameter can be determined independently to high precision, 
from the phase information. 

3 For the instrum ental noise, instead of the approximation of 
Barack & Cutler (2004), we use the more exact sensitivity curve available 
at http://www.srl.caltech.edur shane/sensitivity/ 



0.1 - 



0.01 




M/M 

FIG. 1. — Relative signal-to-noise ratio (SNR) for LISA detections of 
the in-spiral phase of equal-mass supermassive black hole coalescences, as a 
function of total mass, M, and redshift, z- A 1-year observation is assumed 
and the normalization is for M = 2 X 10 6 Mq and z. = 1. 



is much longer than the inverse of the signal frequency, it is 
plausible to assume that the angular position uncertainty de- 
couples from the intricate waveform, and improves linearly 
with the signal amplitude. This assumption is consistent with 
Fig. 1 of Berti. Buonanno. & Will (2005b) which shows that 
the inclusion of spin-orbit and spin-spin terms modifies the 
waveforms but does not alter the angular resolution. In fact, 
the angular resolution is shown to be independent even for 
alternate theories of gravity (e.g. scalar-tensor and massive 
graviton theories). Thus, similar to equation Q, we estimate 
the mass and redshift dependence of the positioning solid- 
angle error as 



tt)(Af,Z): 



SNR(M,z) 
SNR(M ,zo) 



<5fi(M ,z ), 



(3) 



where the -2 exponent assumes that the uncertainty in this 
two-dimensional quantity is the product of independent un- 
certainties in the one-dimensional azimuthal and polar angles. 

A limitation of the above analysis at the high-mass end 
of the range of SMBHs is that these events may not be vis- 
ible for a full year, due to the low-frequency noise wall, be- 
low which LISA looses sensitivity. For instance, for M z > 
9 x 1O 6 M (1.3 x 10 6 M Q ), the low frequency noise wall at 
0.03mHz (O.lmHz) is crossed less than 0.25yr before reach- 
ing the ISCO. At higher masses or redshifts, therefore, the an- 
gular information, which is inscribed in a modulation with a 
lyr period, becomes significantly harder to disentangle from 
other parameters, such as d^, and the errors estimated from 
the SNR alone by equations Q and (j3ji become less accurate. 
In this regime, a better approximation to the scaling of the 
errors is oc (AT /I yr)" 1/2 x SNR' 1 , where AT < 1 yr is the 
time elapsed from the moment the binary appears at the low- 
frequency wall to the moment it reaches the ISCO (H ughes 
2005, private communication; see lHolz & Hughes1d200 5) for 
a more detailed treatment and discussion). 

3. LOCALIZING THE COUNTERPARTS 

We next consider how to use the three-dimensional spatial 
localization of the SMBH merger event by LISA. The solid 
angle error box directly yields the two-dimensional angular 
position error on the sky, in which any EM counterpart will 
be located. However, an additional step is necessary to con- 
vert the luminosity distance, d^, measured by LISA into a red- 
shift, z, which is the relevant observable for the EM coun- 
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terpart. With a particular choice of cosmological parameters 
Pi = (Hq,£Im,Qde,w), we can directly convert a redshift to a 
luminosity distance (see eq. |4] below), and vice-versa. One 
may then envision the following strategy: given the precision 
with which p, are known from other observations, one can 
estimate the redshift, and restrict the search for counterparts 
within the redshift shell corresponding to the L/SA-measured 
luminosity distance, d^. If a counterpart is uniquely identified 
within this redshift shell, and its redshift can be determined 
precisely, then one could hope for an improved measurement 
of the c/l(z) relation, and hence a refined determination of the 
cosmological parameters. 

The first step in this exercise is to determine the expected 
redshift of the source. Apart from the errors on the cosmo- 
logical parameters and on the measurement of d^, the pecu- 
liar velocity of the source relative to the Hubble flow, and its 
magnification due to weak gravitational lensing by inhomo- 
geneities in the distribution of mass along the line of sight, 
introduce two additional sources of errors. In summary, the 
redshift uncertainty will thus include a combination of uncer- 
tainties from (i) the LISA luminosity distance, (ii) the cos- 
mological parameters, (iii) peculiar velocities, and (iv) weak 
g ravitational lens ing magnification. 

iHughes I J2002I) made a simplified estimate of the redshift 
error, without peculiar velocities or gravitational lensing dis- 
tortions, assuming a flat cosmology with a cosmological con- 
stant (assuming VLm = 1 — J7a and w = — 1), and ignoring corre- 
lations with other cosmological parameters. Here, we extend 
that study by using a general form of dark energy (relaxing the 
w prior), by taking into account the various parameter corre- 
lations, and by including errors due to peculiar velocities and 
gravitational lensing. 

To begin, we recall the luminosity distance to a source at a 
fixed comoving coordinate in a smooth Friedmann universe, 

d h (z, Pi ) = (l+z)c T (4) 
Jo H(z',Pi) 

where 

H(z, Pi ) = H ^/n M (l+z)i + n DE (l+z) 3( - 1+K l (5) 

We ignore spatial curvature and set fi^ = 0, in line with pre- 
vious studies and as suggested by recent WMAP data. For a 
source with a small but non-zero radial peculiar velocity, v, 
equation @ is modified, and the luminosity distance is given 
by 

d h {z,pi,v) = d h [z v ,Pi,0], (6) 

where Az = z„— z = (1 +z)v/c is the additional apparent red- 
shift due to the peculiar motion. In an inhomogeneous uni- 
verse, sources along different lines of sight can suffer differ- 
ent amounts of gravitational lensing magnification /i. If [i 
denotes the magnification of the signal power, then the GW 
amplitude, and thus the inferred value of d£\ scales as p, 1 ' 2 . 
For a line of sight that suffers a magnification \i, equation (|6j 
is modified as 

ddz,Pi,v,fi)= d h [z v ,Pi,0,0]. (7) 

In the limit of weak lensing p}l 2 ~ 1 +k, where k <C 1 denotes 
the weak lensing convergence (see below). 
Taking the variation in both sides of equation Q, we obtain 

5d L = -d h K+-—6z+-^-v+) -—dpi, (8) 

az ov ' a pi 
i 



Solving for 5z, taking the square and the expectation value 
of both sides, and using the fact that the LISA measurement, 
the peculiar velocities, and the cosmological uncertainties 
are independent, i.e. (5pt8d{) = (vSdi) = (k<Wl) = { v $Pi) = 
(n5pi) = (kv) = for all i, we find 

(& 2 ) = (^) 2 «<«l.l.sa} + (Sdt, cosm ) + (Sdt >pec ) + (6dl M )) 

(9) 

or equivalently, 

(& 2 > = (&LISA> + (fcLm) + (&pec) + (^wl) , (10) 

where the notation was introduced to distinguish the intrin- 
sic LISA measurement error, (Wl.lisa, from the error resulting 
from cosmological parameters (3d^ cosm ), peculiar velocities 

(Sd^y), and weak lensing magnification (<Wlwi)- We now 
discuss each of these terms, whose forms and magnitudes fol- 
low directly from equation Q. 

3.1. Cosmological uncertainties 
The cosmological term in equation l|9} is 

(wL os J = Ef t ^^^) 5 (id 

~f d Pi OPj 

where d^ and its derivatives are to be evaluated using equa- 
tion @, at the fiducial values of the cosmological parameters 
Pi, and for v = n = 0. In order to place LISA in the context of 
other experiments planned in the next decade, we compute 
('Wlcoshi) from the covariance matrices (SpiSpj) expected 

from two future cosmological probes: Planck 4 (assumed to 
have measured the temperature and polarization anisotropics 
of the cosmic microwave background), and the Large Syn- 
optic Survey Telescope, LSST 5 (assumed to have measured 
the power spectrum and redshift distribution of ~ 100,000 
galaxy clusters). We adopt th e forecasts for the Fisher matri- 
ces for these experiments by Wang et al. (2004). Their analy- 
sis assumes a flat background universe with 6 free parameters 
(Cl DE ,u) M ,w,u) b ,n s ,os) where uj m = Vt M h 2 = (l-il DE )h 2 de- 
fines the Hubble parameter. Note that the luminosity distance 
does not explicitly depend on n s and a%, but they are included 
here, since they couple to the other four parameters as deter- 
mined by Planck and LSST (and hence increase the uncertain- 
ties on these other four parameters). The fiducial parameters 
are (0.7 3,0.14,-1,0.024, 1 0.9), respectively (consistent with 
WMAP; Sperg el et alJ l2003). Since the two observations are 
independent, we simply sum up the two individual Fisher ma- 
trices; the covariance matrix, (SpiSpj), is obtained by taking 
the inverse of the Fisher matrix. In order to substitute in equa- 
tion il Q , it is necessary either to revert to the original cosmo- 
logical parameters (Hq, FIde, ^m, w) in the correlation matrix, 
by performing an orthogonal transformation in the parameter 
space 6 , or to simply write dL(z,pd of equation ® in terms of 
the parameters w m ,w) and evaluate the partial deriva- 

tives in equation il It as a function of these parameters. Fol- 
lowing either approach, we find (^Lcosm) 1 /^l = 1.7 x 10~ 3 
for z = 1. For comparison, we performed the same anal- 
ysis using the Fisher matrices of WMAP (using only tem- 
perature anisotropies) and SDSS (using the power spectrum 

4 See www.rssd.esa.int/index. php?project=PLANCK 

5 See www.lsst.org 

6 The parameters SQ,m and 5Q.de will be fully anti-correlated, because of 
the assumption of flatness. 
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FIG. 2. — Errors on the inferred redshift of an electromagnetic counterpart 
to a LISA coalescence event, for m\ =m.2= 10 6 Mq . The intrinsic LISA error 
on the luminosity distance, d^, is shown as a solid line for the two searings 
5d L /d L oc SNR- 1 ("LISA-a") and Sd L /d L oc AT' 1 / 2 X SNR~ l ("LISA-b"). 
EiTors due to the peculiar velocity of the source (for v = 500 km s _1 ; short- 
dashed line), uncertainties on the background cosmology (long-dashed line), 
and weak lensing magnification (dash-dotted line) are also shown (see text 
for details). 



density perturbations, multiplied by the linear bias of the ha- 
los. Under this assumption, we find that for fixed halo mass, 
the typical peculiar velocity evolves very little from z ~ to 
Z ~ 1 . T his conclusion is con sistent with the semi-analytic 
model of Hamanaet al. (2003), which shows essentially no 
evolution (only a modest decrease from z ~ to z ~ 0.8). 
Therefore, we substitute v ~ 500 km s" 1 in equation H2\ 
at z ~ 1, yielding an error of (Sdl^ c ) l / 2 /d L = 4.1 x 10" 3 . 
The corresponding redshift r.m.s. error contribution (equa- 
tion 1101 ) for v ~ 500 km s" 1 is shown as a function of red- 
shift in Figure[2](short-dashed curve). The figure shows that 
the peculiar velocity error is lower than the typical LISA error 
forz>l. 

Note that the r.m.s. peculiar velocity equals ~ 500 km s -1 
in the most typical cases, but its exact value depends on 
the specific mass of the halo, Mhaio, embedding the quasar 
( Ham ana et al.ll2003h . smaller Mh a i implying lower veloci- 
ties. For a specific source, Mhaio can be estimated directly, 
using the n umber of galaxies that cluster around the identi- 
fied quasar (Kauffmann & Haehnelt 2002). The r.m.s. pe- 
culiar velocity error may then be estimated for this particu- 
lar source, and may be somewhat lower/higher than the value 
shown in Figure^ depending on whether the source is located 
in a galaxy-poor/galaxy-rich environment. 



of the luminous red galaxies, together with the galaxies in 
the main SDSS survey, and following iHu & Haimanl 120031 
for redshift binning, mass limit, and sky coverage). We find 
( 5d t,easJ 1/2 Ml = 1 • 1 x 10" 2 for z = 1 in that case. 

The luminosity distance error corresponds to a redshift er- 
ror according to equation (|9}. Note that equation (II It de- 
pends on the fiducial redshift through the dh derivatives. This 
dependence of the redshift error is shown as a long-dashed 
curve in Figure|2] along with other sources of redshift errors. 
The luminosity distance at z ~ 1000 is measured very accu- 
rately by Planck, and its evolution is essentially unaffected 
by the cosmological parameters down to dark-energy domi- 
nation at z < 2. The figure shows that, as a result, the relative 

cosmological error (fo 2 osm ) 1 ' 2 /z reaches a constant value be- 
yond z > 2. The figure also shows that the cosmology error 
becomes smaller than the typical LISA uncertainty at z > 0.7. 
We find that, even for our best-case LISA events, the cosmol- 
ogy error becomes sub-dominant at z > 1 . 

3.2. Peculiar velocities 

Equations (|6j, (|8j, and (|9j yield, to first non-vanishing order 
in v/c, the uncertainty due to peculiar velocities, as a simple 
function of the r.m.s. peculiar velocity of the GW sources, 
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We are assuming here that GW events correspond to lumi- 
nous quasars, and unfortunately, the r.m.s. peculiar velocity 
of high-redshift quasars is not known empirically. We there- 
fore employ theoretical predictions for the peculiar veloci- 
ties of (quasar-host) galaxies within d ark matter haloes. Ac - 
cording to numerical simulations (e.g. Cen & Ostriker 2000), 
typical values are ~ 500 km s -1 , with a tail extending to ~ 
1000 km s . As an approximation, we assume that the scal- 
ing with redshift follows the linear growth in the amplitude of 



3.3. Weak Gravitational Lensing 

The gravitational lensing term in equation (|9jl, in the weak- 
lensing limit, is given by 



(SdlJ/dl- 



(13) 



where k denotes the r.m.s. effective convergence (e.g., 
iWhite & HiJ l2000). While the mean magnification is well 
approximated by (k) = or (/i) = 1, the magnification dis- 
tributi on has a su bstantial w idth. T he full distribution is given 
by Wang, Holz, & Munshi (2002), and its variance reaches 
~ 12% for sou rces at 7 = 2 llDalal et alJl2003l) . Here we use 
equation (6) in Whi te & Hul J2000I) to compute the variance 
in effective convergence, (re 2 ), for point sources. This quan- 
tity is given by an integral over the matter power spectrum, 
and receives a contribution from small, non-linear sc ales. We 
employ the HALOFIT routine o flSmith et. al. M2003I) . and set 
the input cosmological parameters according to our conven- 
tion (see §Q. This routine encodes an accurate fitting formula 
for the matter power spectrum extending into the nonlinear 
regime. For z = 1, we find (k 2 ) 1 / 2 = 3.1%. 

Weak lensing errors can, in princ iple, be improved by "cor- 
rective lenses' i.e. background galaxy 
shear maps, and using the cross-correlation between the mag- 
nification of a point source and the sh ear m ap smoothed on 
larger-scales. However, iDalal et al] (|2003) found that the 
magnification errors can be improved by only a small amount, 
less then 20% relative to the uncorrected errors for a source at 
z = 2. 

A different approach for reducing weak lensing magnifi- 
cation uncertainty would be to directly measure the inhomo- 
geneities in the mass distribution along the line of sight. If 
this distribution could be directly probed down to a scale of 
k m \ n , then the contributions to K from all scales down to k m \ a 
could be directly subtracted from the uncertainty on the mag- 
nification. If the counterpart is indeed a quasar, then the line- 
of-sight density distribution could, in principle, be probed by 
studying its Lyman a absorption spectrum, as well as deep 
surveys of galaxies and clusters in the foreground and near the 
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FIG. 3. — The relative improvement in the weak lensing magnification- 
induced redshift uncertainty to a background source. We show the fractional 
improvement that can be achieved by perfectly measuring density inhomo- 
geneities down to a fixed scale. The .v-axis shows the minimum wave number 
that is unmapped (larger scale fluctuations are assumed to be perfectly known 
along the line of sight to the source). Various curves correspond to sources at 
Z= 1, 2, 3, and 4 (from bottom to top). 



line of sight. At low redshifts (z ~ 0.5) which contribute sig- 
nifican tly to the lensing magnification, the X-ray abs orption 
forest ( Fang & Canizares 2000; Perna & Loeb 1998) could 
provide additional information on the density fluctuations. 

We leave a detailed assessment of the amount of correc- 
tion that could be feasible to future investigations. However, 
to estimate the "target" scale at which the fluctuations would 
need to be measured in order to be useful as a lensing cor- 
rection, we make the (unrealistically optimistic) simplifying 
assumption that the matter fluctuations above a certain length 
scale, and the corresponding contribution to lensing magnifi- 
cation, have been perfectly determined. The fluctuations on 
smaller scales, with k > & m ; n , then determine the only remain- 
ing weak-lensing uncertainty. Therefore, we truncate the in te- 
gral over the wavenumber (see eq. [1] in Dalai et al. 2003) at 
^min- The resulting fractional improvement in (k 2 ) 1 / 2 (which 
equals Sdh,wi/d]_) is plotted in Figure [3] The improvement is 
about 8%, 20%, and 40% for k min =1,3, and lOMpc" 1 , respec- 
tively - i.e., about half of the weak lensing uncertainty is from 
length-scales < 27r/£: m j n ~ 0.6 Mpc. The Lyman a transmis- 
sion spectra of SDSS quasars have been used to determine the 
power spectrum down to scales of k ~ 3 Mpc -1 , implying that 
the spectral resolution would need to be improved by a factor 
of ~ 3 to allow significant improvements. 

The uncorrected weak-lensing redshift uncertainty 5z w i = 
(dd L /dz)~ l {Sd^ wl ) l l 2 is shown in Figure [5] The plot clearly 
shows that weak lensing errors typically exceed 1% and dom- 
inate the other redshift errors for z > 0.5. 

3.4. The Size and Orientation of the Error Volume 
In Figure |2j we compare the contributions from four 
different sources of redshift errors: (&ljsa) > (^ z pec) ' 
(&co.sm> 1/2 - and (&wi>' /2 for SMBH masses m x = m 2 = 
10 6 M Q . Weak lensing errors are dominant for 5 > z > 0.5. 
The LISA redshift error is shown for both scalings Sd^/d^ oc 
SNR~ l (labeled "ZJSA-a") and 6d L /d L oc (AT/1 yr)" 1 / 2 x 
SNR' 1 (labelled "LISA-b"). These curves separate for M z > 
3.92 x 10 6 M Q , when the maximum observation time AT de- 



creases under 1 yr. The other three sources of errors are within 
a factor of two of one another at z = 1 , for typical LISA sources 
and peculiar velocities. 

The total redshift uncertainty then follows from summing 
the uncertainties in quadrature (eq. 1101 ). The result for 
„H = m2 = 10 6 M Q and z = 1 is Sz = 2.57%, 2.59%, and 3.03% 
in the cases we labeled as best, typical, and worst, respec- 
tively. Note that we have not explicitly added an error due to 
the instrumental resolution of a spectrograph, because this is 
often much better than this value (e.g. for SDSS 7 , the redshift 
resolution is between 10~ 3 and 10~ 4 ). 

We next use these redshift uncertainties to derive the three- 
dimensional error volume, SV Xot , in which the EM counterparts 
need to be identified. The comoving volume corresponding 
to the above redshift uncertainties, combined with the solid- 
angle uncertainty (SCl 2 ) 1 / 2 , is 

^=S (&2>1/2( ^ 2>1/2 - (14) 

Substituting the solid-angle errors from Table [I] we find 
5V M = 2.1 x 10 3 , 6.3 x 10 4 , and 7.4 x 10 5 Mpc 3 for the same 
three cases. 

Equation J14i should be regarded as a simple estimate of 
the volume that needs to be searched. The 2D uncertainty 
((Jfi 2 ) 1 / 2 was obtained by a Fisher analysis, and represents 
the solid area of a 2D error ellipse. This equation then gives 
the volume of a cylinder with a height of (fo 2 ) 1 / 2 . The ar- 
guably more appropriate volume of an ellipsoid whose third 
semi-axis is (Sz 2 ) 1/2 would be a factor of 4/3 larger. On the 
other hand, equation dl4> also naively assumes that Sz and 5 ft 
are uncorrected - i.e., it describes an error ellipsoid whose 
z axis is oriented along the line of sight. The angular po- 
sition is estimated by LISA from the GW signal alone, and 
is therefore indeed uncorrelated with the radial uncertainties 
due to cosmology, lensing, or peculiar velocities. However, 
the luminosity distance estimate from LISA itself is strongly 
correlated with its angular positioning (Hushes 2002), which 
results in an error ellipse that is 'tilted' relative to the line of 
sight, and has a smaller overall volume than the simple or- 
thogonal product in equation (II 41 would imply. 

We have utilized the correlation matrices for the LISA dis- 
tance a nd angle estimates given in Table 1 and 2 of Hughes 
J2002I) for mi=m 2 = 1O 5 M at z = 1 and m, = m 2 = 10 4 M Q 
at z = 7, to estimate the reduction in the total error volume 
due to these correlations. Note that this analysis applies only 
to the LISA un certainties. Among the 1 1 free parameters in 
Hughes I (12001 . the 3D error volume is determined by the pa- 
rameters related to the spherical coordinates of the sources: 
lnc/L, A*5 = cos 9s, and 4>s- The error volume SV to t quoted in 
equation dl4> above corresponds to (3/4th of) the volume of 
an ellipsoid whose semi-axes are the marginalized errors in 
spherical coordinates; the true error volume SV e \\ is that of the 
ellipsoid described by the full covariance matrix. The ratio 
SV e \\/SV tot depends on the actual position angle (0s, <As) of the 
source; averaging over all angles, we find (<5V e u/<5Vtot) = 0.31 
(0.20) for masses m x =m 2 = 1O 5 M (10 4 ) at redshift z = 1 (7). 

If LISA errors dominated the total redshift uncertainty 
(fo 2 ) 1 / 2 , this implies that the correlations could reduce the 
mean number of counterparts by a factor of 3^1. However, 
as discussed above, the total uncertainty is likely going to be 
dominated by weak lensing errors; hence the inclusion of the 

7 http://cas.sdss.org/dr3/en/tools/search/sql.asp 
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correlations would reduce the final error volume only by a 
small factor (~ 15%). We make no use of this reduction in 
the results we quote below. 

4. QUASAR COUNTERPARTS 

To estimate the typical number of quasar counterparts to 
a specific SMBH merger event, we combine the size, 5V co , 
of the comoving LISA error box with the space density of 
quasars, by integrating over the quasar luminosity function, 

6Lcj)(L,z), (15) 

where b accounts for the bias due to the clustering of quasars, 
and L m ; n and L max are the minimum and maximum quasar lu- 
minosities which could be associated to the specific SMBH 
merger event. We use = O.lLEdd and L max = 2L^&&, where 
^Edd denotes the Eddington luminosity for the total BH mass. 
Motivations behind this particular near-Eddington choice are 
further discussed below. 

4. 1 . Luminosity Function of Quasars 

We adopt the standard empirical double power-law fit to 
the quasar luminosity function of the combined quasar sam- 
ples from the Two-Degree Field (2dF) and Sloan Digital 
Sky Survey, with pure luminosity evolution valid for z < 2.1 
( Rich ards et a1~l 120051 ICroom et all 120051) . Unfortunately, 
these surveys extend only to relatively bright magnitudes (B < 
21), corresponding to the Eddington luminosity of BHs with 
mass M > 3 x 10 7 M Q for z = 1, which is above LISA'S op- 
timal mass-range. In order to estimate the number of quasar 
counterparts for lower BH masses, we extrapolate the lumi- 
nosity function to the faint end. 

A limitation of this luminosity function is the si mple 
quadratic fitting formula for the evolution dRichards et al. I 
120051 Croom et al. 2005), which is only valid for re dshifts 
z < 2. 1. To avoid these difficulties Madau. Iiaardt. & Reesl 
( 1999) (hereafter MHR) proposed a more complicated em- 
pirical fitting formula with three adjustable parameters z s , C> 
and £. These parameters were estimated using high-z quasars, 
and were improved to include the high-redshift SDSS sam- 
ple and weak-lensing effects ( Wvit he & Loebl2002l) . We ob- 
tain a luminosity function that is more precise at lower red- 
shifts and is conco rdant with the re dshift scaling by fitting the 
MHR model to the Richa rds' et al. I (120051) luminosity function 
in the interval .5 < z < 2.1 and keepi ng the high-redshift 
asymptote at the I Wvithe & Loebl J2002) value. The result is 
L* = 5.06x 10 10 L Q ,z s = 1.66, C = 2.6, and £ = 2.8. Other pa- 
rameters in the luminosity f unction that are i ndependent of the 
evolution are adopted from Richar ds et al. 1 (12005): (3\ = 1.45, 
f\ = 3.31, and $* = 1.99 x 10- 6 Mpc" 3 . 

4.2. Clustering of Quasars 

The bias, b, in equation dl5> describes the enhancement in 
the number of quasars around a specific quasar being the po- 
tential counterpart, relative to the value for a homogeneous 
distribution . The clustering depends only weakly on quasar 
luminosity (iLidz et al.l2005t lAderbei-ger & Steide ll2005l) . We 
use the observe d autocorrelation function of quasars from 
the 2dF survey ( ICroom et alJ 120051) . £(s) = (s/so) 7 for s > 
0.1/i _1 Mpc (assuming no quasars with a smaller separation), 
where s Q = 5.48/r'Mpc (5.55) and 7 = 1.20 (1.63) for s < 



TABLE 2 
Survey characteristics 



Survey 




Limiting mag 


Sky 


cov. 


Observing 




(z=l)M Q 




de. 


g 2 






2dF a (2QZ.6QZ) 


3 x 10 7 


bj < 20.85 






750 


1997-2002 


SDSS b (LRG) 


2x 10 8 


(< 19.1 


7 


X 


10 3 


1998-2005 


Deep2 c 


5 x 10 5 


R < 24.5 






4 


2002-04 


AGES d 


1 x 10 7 


i?<21.5 






9 


2004-06 


DES C 


7 x 10 5 


AB < 24.7 


5 


x 


10 3 


2009-13 


LSST f 


2 x 10 5 


AB < 26.5 


1.8 


X 


10 4 


2012 



a Two-Degree Field, see http:/www.2dfquasar.org/ 
b Sloan Digital Sky Survey, see http:Avww.sdss.org/ 
c See http:/deep. berkeley.edu/ 

d AGN and Galaxy Evolution Survey, covers radio, IR, optical, and X-ray 
bands, see http:/cmb.as. arizona.edu/~eisenste/AGES/ 

c Dark Energy Survey, see http:/cosmology.astro. uiuc.edu/DES/ and 
http :/decam . gnal . gov/ 

f Large Synoptic Survey Telescope, see http:/www.lsst.org/ 



25Mpc (>25Mpc), respectively. The bias is given by the aver- 
age value of (1 +£(?)) over the comoving error box. Assum- 
ing that the error box is a cylinder, with height Sy = c5z/H(z) 
(which is the comoving distance along the line of sight corre- 
sponding to the redshift error) and radius 5r = \/SV C0 /{Tr5y) 
(corresponding to the angular uncertainty of LISA): 

f* r 2nrdrf^dy(l + aVy^)) 

b= fx T' ' ' > 

J x 2nrdrf y dy 

We find b = 1.50, 1.23, and 1.07, for our best, typical, and 
worst cases, respectively. The corresponding comoving cylin- 
der heights are 5y = 62.5, 62.9, and 73.8 Mpc, and the cylinder 
radii are Sr = 3.3, 17.8 and 56.4 Mpc. 

5. RESULTS 

The main results of this paper are presented in Figure |4] 
which shows (N), the average number of counterparts within 
the 3D LISA error volume, for various total masses M and 
redshifts z. Recall that we have assumed that the GW event is 
always accompanied by quasar activity. According to our def- 
inition, (N) then corresponds to the mean number of quasars 
that would be found in LISA'S error box, in addition to the 
quasar actually associated with the GW source. A straight- 
forward identification of a unique counterpart therefore re- 
quires that there be no additional quasars in the error vol- 
ume and that the EM observation sensitivity goes below the 
actual counterpart luminosity (so that the presence of fainter 
quasars can be ruled out). A simple criterion for a reason- 
able chance not to have any additional counterparts candi- 
dates is (N) < l. 8 We find that this simple condition is sat- 
isfied in the case of a "typical" event at z = 1 with total masses 
~ 4 x 10 s M© or — 8 x 10 6 M Q . At higher redshifts, the av- 
erage number of potential counterparts will be much larger, 
due mostly to the increasing weak lensing errors. At z > 3, 
even the best case events will typically have at least one ad- 
ditional quasar in their error box. On the other hand, the 
increase from z = 3 to z = 5 in the number of quasars lo- 
cated in the error-box is partly mitigated by the drop in the 

8 One could explicitly consider the probability distribution for N. For ex- 
ample, for a Poisson distribution, the probability for N = would be 50/90 
percent for (N) =0.7/0.1. 
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abundance of quasars at z > 3. The three panels of Fig- 
ure 0] display results for various presumptions of uncertain- 
ties. The top panel uses raw, uncorrected weak-lensing er- 
rors and the counterpart luminosity is allowed to vary in a 
broad range 0.1 < L/Lecu < 2. The middle panel accounts 
for a 20% reduction of weak-lensing errors, and the lumi- 
nosity is restricted to 0.7 < L/L E dd < 1.3. The bottom panel 
shows results for the more conservative scaling of LISA er- 
rors, 6d L /d h oc Ar'/ 2 x SNR~ l and dQ, oc AT~ l x SNR' 2 if 
AT < lyr. In all cases at z ~ 1, a unique counterpart may be 
identifiable. 

One may also interpret the number of counterparts we com- 
pute, (N), together with the same simple criterion (TV) < 1, in 
a somewhat different way: as a means to test our hypothesis 
that LISA events are accompanied by bright quasar activity. 
If our hypothesis is incorrect, then the number of quasars in 
LISA'S error volume should be drawn from the random dis- 
tribution of quasars on the sky, unrelated to the LISA event 
(excluding the correction due to correlations that we have in- 
cluded in our analysis; although we found this correction to 
be relatively insignificant). In many configurations, we find 
(N) <C 1, implying a significant probability that no bright 
quasars would be found in LISA'S error volume in these cases. 
If several LISA events are indeed found with no such quasar 
counterparts, it would, by itself, be an important new con- 
straint on the process of binary black hole coalescence. 

The identification of counterparts could be aided by a com- 
bined analysis of several GW events. Every successful iden- 
tification yields a very precise direct measurement for the 
L/Lem Eddington ratio. Once a statistically significant set 
of Eddington ratios is acquired, the empirical distribution can 
be mapped. It is yet unclear whether SMBH mergers are ex- 
pected to have high-luminosity quasar counterparts. Then, in 
case L/Lem is in a narrow range, this information can be used 
to greatly constrain the a priori assumptions on the counter- 
part's luminosity. As an example, in Figure |4] (middle panel) 
we consider L/Lem = 1 ±0.3, and find (TV) to decrease well 
under 1 in the typical case. 

If cosmological uncertainties were to dominate the error 
budget on a counterpart's redshift, a combined analysis could 
further improve the robustness of the identification. Indeed, 
even if each GW event has, by itself, several possible coun- 
terparts, each of these counterparts would require a different 
set of cosmological parameters. As a result, there should be, 
in general, only a single set of cosmological parameters 9 that 
gives a con sistent set of counterparts to all of the GW events 
(Schutz 1986). The counterpart candidates contradicting this 
concordant set can be discarded. Unfortunately, the error bud- 
get is likely to be dominated by the lensing magnification un- 
certainty. In this case, having multiple events is going to be 
helpful only if a sufficient number (^ 100) of events are de- 
tected to map out the full distribution of magnifications (such 
a large number of events is not expected; see below). 

6. DISCUSSION 

6.1. Search Strategy 

Different strategies for the search of an electromagnetic 
counterpart to a LISA event can be envisioned. The simplest 
one would be to search, in an existing survey catalog, for can- 
didates located in the 3D LISA error box. We list the char- 

9 This set, of course, will suffer from the usual degeneracy along the sur- 
face of constant (2l(z, Q\,Q m ,h,w) around the fiducial cosmological param- 
eters. 
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Fig. 4. — The average number of quasars in the three-dimensional LISA 
eiTor volume. The thin dash-dotted, thick solid, and thin dashed curves trace 
the "best", "typical", and "worst" cases, respectively. The low-frequency 
noise wall for LISA is assumed to be at f m \ n = 0.03mHz. Top: Using raw 
data without any weak-lensing corrections and the counterpart's luminosity 
is assumed to be 0.1 < L/Lridd < 2. Middle: The weak-lensing errors are 
corrected by 20%, and the luminosity is assumed to be 0.7 < L/^Edd < 1-3. 
Bottom: Same as the middle panel, assuming that LISA uncertainties scale 
as Sd L /d L oc AT' 1 / 2 X SNIT 1 and dtt oc AT' 1 X SNIT 2 . In most cases at 
Z < 1, a unique quasar counterpart may be identifiable. 
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acteristics of several deep present and near-future surveys in 
Table|21 These surveys cover various bands between 400 and 
llOOnm. In each case, we calculate the minimum mass of 
the quasar corresponding to the magnitude limit of the survey, 
provided that the quasar is shining at the Eddingt on luminos- 
ity. Fo r the quas ar spectrum w e use Madau. liaardt. & Reesl 
(1999) and follow M adau I i ll 9951) to account for the Lyman se- 
ries line blanketing and photoelectric HI absorption to obtain 
the apparent luminosity. The portion of the sky covered and 
scheduled operation times are also given. The most promising 
(deepest and largest sky coverage) among the planned surveys 
is that from the LSST, scheduled for the same time as LISA. 

A second strategy would be to design a survey specifically 
aimed at identifying the counterparts of well-localized LISA 
events, by observing in detail the 2D angular error box pro- 
vided by LISA and looking for the expected counterparts. The 
time domain may come to play an important role in the strat- 
egy of this type of surveys. Indeed, for some of the events, 
one may have, from the LISA data-stream, a reasonably good 
idea of the location on the sky of the SMBH bina ry prior to 
coalescence. For instance, Holz & Hushes (2005) estimate 
that this knowledge may be available with reasonable accu- 
racy about a day in advance. One would then be able to 
monitor any unusual photometric variability associated with 
the violent SMBH merger in this area somewhat before, dur- 
ing and after the coalescence. Surveying the area long (i.e. 
months or years) after the coalescence may also prove useful 
in discovering the counterpart and/or monitoring the viscous 
evolution of any gas surro unding the SMBH merger remnant 
(Milosavlievic & Phinnev 2005). The time domain may thus 
greatly facilitate the identification of a unique electromagnetic 
counterpart to some LISA events, even in cases when many 
counterpart candidates otherwise exist in the 3D error box. 

6.2. Uncertainties in the Analysis 

Throughout the analysis of LISA uncertainties, we focused 
on equal mass binaries. With this simplifying assumption, 
it was possible to derive approximate scalings as a function 
of tota l mass and relate them to the calculations of IVecchiol 
(2004) for m\ = '«2 = 1O 6 M0. We also find good agreement 
with the recent calculations of Holz & Hughes (2005) for 
a variety of equal-mass combinations. For unequal masses, 
the GW signal-power depends to leading order on the simple 
combination M. = /i^M 2 / 5 , allowing a straightforward ex- 
tension of our scaling arguments to more general mass ratios. 

As mentioned above, our results depend on the LISA sensi- 
tivity curve at low frequencies. In particular, the relationship 
between the arrival frequency f a {M z , AT), the final frequency 
/i sco (M : ), and the LISA minimum frequency noise wall / m ; n 
determines the maximum possible observation time for the 
in-spiral. The exact value of fmm is currently assumed to lie 
between 0.1 and 0.03mHz. Figure|4]assumes fmm = 0.03mHz. 
In this case, at z = 1, a AT = lyr observation of the in-spiral 
phase is possible for M < 2 x 10 6 M Q but for larger masses, 
the maximum possible observation time becomes shorter than 
a year. We have also computed the number of quasars in the 
3D LISA error box for the higher value of / m ; n = 0. lmHz. We 
find a significant increase in this number at the high-mass 
end (i.e. forM ~ 10 7 M Q ), but the results are essentially un- 
affected for lower mass SMBHs. 

An important assumption in our analysis is the near- 
Eddington luminosity of the quasar counterparts associated 
with LISA events. While this assumption is difficult to justify 
from first principles, it is the luminosity expected if SMBH 



accretion occurs in a regime which is not supply-limited. Ob- 
servationally, the Eddington ratio of quasars, L/L^, can be 
inferred, leading to values fro m 0.1 to > 1, with higher val- 
ues a t large redshifts (e. g. iKaspi et ail 120001 IVestereaardl 
2004; Woo & Urrvl 12002). This ratio has been determined 
for a handful of lower mass SMBHs in AGNs and these 
sources have been fo und to cluster aro und L/Lfm ~ 1 as well 
( iGreene & Hol2004l). Rec ent work byHollmeier eTal] (120051) 
and Hopkins et al. (2005) both suggest L/Lndd is typically 
around 1/3. 

The restricted luminosity range 0.1 < L/Lgdd < 2 assumed 
in our analysis effectively serves as a 4th ("brightness") di- 
mension, complementing the three-dimensional geometrical 
error volume provided by LISA. According to Tabled for a 
given event, the mass estimate provided by LISA is typically 
very accurate. Therefore, it is through the range of acceptable 
Eddington ratios that the integration bounds in equation Jl 51 
change. Since the luminosity function decreases rapidly with 
L, the integral is dominated by the lower bound and is largely 
independent of L mix . We chose L max = 2L£dd for concrete - 
ness. When modifying the lower bound from L m j n = 0. ILecM to 
Lmin = 1 OLEdd or L m i n = O.OlLEdd, for instance, we find that the 
number of quasars present in the 3D LISA error box changes 
by a factor of 0.1 and 5, respectively. 10 Obviously, our need 
to extrapolate the quasar luminosity function below the min- 
imum value constrained by the observations is an additional 
important source of uncertainty in our analysis, but one that 
will be addressed well by future generations of surveys (Ta- 
bleEJi. 

In a recent study, Hennawi et alJ d2005l) have found a large 
number of small-separation quasars, implying an order-of- 
magnitude increase in the quasar auto-correlation on scales 
< 200/T 1 kpc. However, the dimensions of the error volume 
as we estimated it is generally much larger (~ 3.3 Mpc even 
in our best case). As a result, the average number of coun- 
terparts is changed by a negligible factor by this small-scale 
auto-correlation: 8.3 x 10" 3 , 3.7 x 10" 4 , and 3.3 x 10" 5 in the 
best, typical, and worst case LISA errors. Increasing the cor- 
relation by an order of magnitude out to a scale of 500/i _1 kpc 
would still cause at most a 6% increase in the m ean number 
of quasars in the best case. Hennawi et al. (2005) do not find 
an increase in the correlation beyond ~ 500/i _1 kpc scales. 
We note that our inability to identify the correct counterpart 
from among any rare ultra-small-separation candidates would 
not degrade cosmological parameter estimations significantly. 
One may also argue that such close-separation binary quasars 
would better represent a 'precursor' stage in the evolution of 
the two black holes towards an eventual coalescence. If so, 
one would not expect them to be associated with GW events, 
and they may in fact be anti-correlated with such events. 

Throughout this paper, we have focused on optical quasars 
as plausible counterparts to GW events. It would be inter- 
esting to repeat our analysis using other types of "electro- 
magnetic objects" that may be associated with SMBH co- 
alescences. For example, even if gas accretion leads to 
prodigious energy output in radiation, the optical light of 
the quasar may be obscured by the intervening gas and dust 
near the galactic nucleus. In these cases, the GW events 
may be more commonly assoc iated with X-ray quasars (e.g. 
Milosavlievic & Phinnev 200 3), or with ultraluminous in- 
frared galaxies (e.g. Thompso n et alJ2005l) . As these sources 

10 The change is explained by the asymptotic form of the luminosity func- 
tion <E>(L) <x ZT 1 ' 45 for low L. 
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are also relatively rare, unique counterparts may be identifi- 
able if such objects typically accompany GW events. 

6.3. Implications: Black Hole Astrophysics 

A successful identification of a quasar counterpart to a LISA 
event would provide powerful diagnostics on the physics of 
SMBH gaseous accretion and the associated radiation. The 
masses and spins of the two BHs before merger can be directly 
determined from the GW signal, from which the mass and 
spin of the remnant BH follows at some precision. In some 
cases (e.g. Hughes & Menou 2005), the remnant BH could 
also be observed by LISA during the post-merger ring-down 
phase, which would constrain its mass and spin directly. The 
orientation of the orbital plane of the BH binary before merger 
would be measured as well. All of these parameters are gen- 
erally unknown for quasars detected only via traditional elec- 
tromagnetic techniques. 

The observation and monitoring of quasar counterparts to 
LISA events may thus offer us some of the best laboratories 
for the study of AGN physics. First, the Eddington ratio can 
be measured to high accuracy (limited only by photometric er- 
rors and bolometric corrections), since LISA estimates for the 
BH masses are extremely precise by astronomical standards 
(see Tab.Q. Second, if the quasar accretes at or near the Ed- 
dington limit, given its Eddington ratio, one may be able to set 
a useful lower limit on the radiative efficiency of its accretion 
flow. If it is in excess of the canonical 10% value, it will pro- 
vide an interesting empirical test of the physics of accretion 
onto a spinning BH, with a spin directly constrained by the 
GW measurement. Third, the counterpart could be monitored 
for years following the merger to follow the viscous evolu- 
tion of the gaseous disk and thus clarify its role in the SMBH 
coalescence process. Fourth, it is expected that the gas disk 
will be forced in the orbital plane of the pre-m erger bin ary by 
the Bardeen-Peterson effect (Milosavlievic & Phinnev 2005). 
Knowing the disk orientation could thus offer tests of the ge- 
ometry of quasar emission and obscuration (even after the 
merger, given the expected spin of the remnant). It may also 
be possible to further develop diagnostics related to the geom- 
etry of a jet, if present. 

6.4. Implications: Cosmology 

The successful identification of an EM counterpart to a 
GW event could, in principle, open the way to use them 
as "standard sirens" to probe the background cosmology 
(Schutz 1986), analogously to the la SNe standard candles 
I Holz & Hughesl l2005l) . The precision on the cosmological 
model can, however, be improved only if the diXz) function is 
determined to a higher accuracy than it can be already guessed 
from other data that exists when LISA is operational. As an 
example, we have assumed here to have available the uncer- 
tainty from the combined datasets from two future projects, 
Planck and LSST We have found that the major obstacle 
against a dramatic improvement on cosmology is the gravi- 
tational lensing of intervening matter along the line of sight 
to a LISA source. In the weak-lensing limit, the r.m.s. mag- 
nification of a source at z = 2 is 12%, lead ing to a luminos- 
ity distance error of 6%. iDalal et all |2003) have shown that 
galaxy shear maps can be used to correct weak-lensing distor- 
tion, but only a 20% relative improvement can be achieved, so 
that (Ml/^l = 5%. We suggest alternatively that correcting for 



the contribution of the known distribution of intervening mat- 
ter might improve weak lensing uncertainties by another 20%. 
Furthermore, the weak lensing uncertainty can be overcome if 
a large sample of sources is available when fitting d^iz), map- 
ping out the magnification distribution. Using a number K of 
merger events along with uniquely identified counterparts, the 
lensing error reduces approximately as 1/ \/~K (although the 
actual improvement will be less pro nounced, due to a non- 
Gaus s ian tail of high magnif ications; IWang. Holz. & M unshi 
120021 [Holz & Linder1 l200l . To reach the Planck + LSST 
level of Sd^/d^ = 10~ 3 , we find that K > 100 events would be 
required. 

Is such a large number of merger events, with uniquely 
determined electromagnetic counterparts, expected from 
the LISA data-stream? Monte-Carlo simulations of 
SMBH merger trees generally indicate LISA event rates 
from ~ 20-0.5yr~ 1 for masses M B h < 10 7 at z ~ 1 
JMenou. Haiman' & Narayanan! 1200 It iSesana et ai~l 12004 
lEnokietalJl2004l) . Detecting a total of K = 100 would be 
barely possible during a 3-year LISA mission lifetime, only 
allowing a marginal test of the concordance of cosmolog- 
ical parameters with standard sirens. In addition, most of 
these events may be expe cted to involve SMBHs at th e low- 
massend(i.e. < 10 5 M ? : lMenoull2003USesana et aT~ll2004l) . 
which are not ideal for unique counterpart identifications (see 
Fig-IU- However, large uncertainti es remain on the e xpected 
event rate. For example, Isla m. Tavlor. & Silkl (|2004) predict 
much larger rates, which could open up the possibility of a 
statistical analysis, folding in the expected weak lensing mag- 
nification distribution. 

7. CONCLUSIONS 

In this paper, we have considered the possibility that 
SMBH-SMBH mergers, detected as gravitational wave 
sources by LISA, are accompanied by gas accretion and 
quasar activity with a luminosity approaching the Eddington 
limit. Under this assumption, we have computed the num- 
ber of quasar counterparts that would be found in the three- 
dimensional error volume provided by LISA for a given GW 
event. We found that weak lensing errors exceed other sources 
of uncertainties on the inferred redshift of the electromagnetic 
counterpart and increase the effective error volume by nearly 
an order of magnitude. Nevertheless, we found that for merg- 
ers between ~(4x 10 5 - 1O 7 )M SMBHs at z ~ 1, the error 
box may contain a single quasar with a B-band luminosity 
Lg ~ (10 10 - 10 n )LQ. This would make the identification of 
unique electromagnetic counterparts feasible, allowing pre- 
cise determinations of the Eddington ratio of distant accret- 
ing SMBHs, and providing an alternative method to constrain 
cosmological parameters. 
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